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ABSTRACT 

We use data from the Pan-STARRSl survey to present a panoramic view of the Sagittarius tidal 
stream in the southern Galactic hemisphere. As a result of the extensive sky coverage of Pan- 
STARRSl, the southern stream is visible along more than 60° of its orbit, nearly double the length 
seen by the SDSS. The recently discovered southern bifurcation of the stream is also apparent, with 
the fainter branch of the stream visible over at least 30°. Using a combination of fitting both the 
main sequence turn-off and the red clump, we measure the distance to both arms of the stream in the 
south. W e find that the distances to the bright arm of the stream agree very well with the N-body 
models of lLaw fc Maiewski (|2010T l . We also find that the faint arm lies ~5 kpc closer to the Sun than 
the bright arm, similar to the behavior seen in the northern hemisphere. 

Subject headings: Galaxy: halo - Galaxy: structure - Local Group 



1. INTRODUCTION 

The past decade has seen a tremendous growth in 
the a mount of known tidal substructure in the Galactic 
naio (Iba ta et a D200lUNewberg et alJ2002HYanny et all 
2001 IMartin et all l2004MGrillmair fc Dionatosl 120061: 
Grillmairl 120061: IBelokurov et all I2007bllal : Uuric et all 
20081) . One of the most prominent and well- 



studied tidal streams in the halo is the Sagittarius 
stream. A wide variety of tracers have been used 
to study the strea m, including main seq uence turn- 
off (MSTO) stars (|Belokurov et all I2007al), blue hor- 
izontal branch (BHB) stars dNiederste-Ostholt et al 
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2010 1 : iRuhland et all 120 111) A st ars (jYannv et al' 



20001: iMartmez-Delgado etail 12001 |200"1, M giants 



Maiewski et all I2003D, R R Lyrae (jVivas et al 



giants 
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2005i I Watkins et all 120091) and red clum p (RC) stars 



(jBellazzini et alll200a ICorrenti et al.H2010f ). These stud- 
ies have all revealed that the Sgr stream is exceed- 
ingly complex and highly structured. The northern 
"bifurcation" , where the stream appears to split into 
two parallel streams, remains unexplained. Multiple 
possible causes of the bifurcation have been proposed 
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(e.g.. IFellhauer et all 120061: IPenarrubia et all 120101 ) but 

no explanation so fa r has proved satisfactory (e.g., 
IPenarrubia et al.ll201lD. 

Recent work by IKoposov et~aT1 (|20ll) using the SDSS 
has shown that the Sgr stream is more complex than 
previously believed, with a bifurcation in the stream in 
the southern Galactic hemisphere as well as in the north. 
This has only added to the challenge of reconciling the 
diverse properties of the stream. It also changes the 
range of possible explanations, e.g., it is possible that 
each of the two northern streams comes from a differ- 
ent progenitor. As a result of all these features of the 
stream itself, along with uncerta inties about the gravita- 
tiona l potential of the Gal a xy (iHelmil 120041 lLaw et all 
[2005t I Johnston et all [20051 IFellhauer etldT 120061 ). the 
Sgr stream has proven to be a serious challenge to 
any model that attemp ts to describe it. The model of 
lLaw fc Maiews ki (2010) has been able to generally re- 
produce many of the features of the stream (assuming a 
triaxial Galactic potential), but the bifurcations are not 
yet accounted for. 

These modeling challenges underscore the importance 
of obtaining a complete set of observations of the stream, 
including data from both the northern and southern 
Galactic hemispheres. In this work we focus on the 
southern component of the stream, using data from the 
Pan-STARRSl project. Pan-STARRS's extensive sky 
coverage allows for a comprehensive view over approxi- 
mately 60° along the southern stream, yielding a picture 
in the south comparable to that available in the north. 

In this work we extend the coverage of the southern 
Sgr stream seen bv lKoposov eTaLl mm in SDSS, nearly 
doubling the length of the stream observed. We measure 
distances along the entire observed southern stream, us- 
ing red clump (RC) stars and a direct calibration of the 
RC absolute magnitude from the Sgr dwarf itself. The 
red clump is an excellent distance indicator, since it oc- 
cupies a narrow range of absolute magnitudes and is only 
weakly dependent on the parameters of the stellar pop- 
ulation. However, the red clump is a relatively weak 



2 



Slater et al. 



feature, and in faint systems it can be difficult to unam- 
biguously identify without additional information. This 
is particularly true at low galactic latitudes where the 
large number of disk stars increases the difficulties of 
identification. To mitigate this problem, we take advan- 
tage of the fact that the main sequence turn-off (MSTO) 
of the Sgr stream is much more well populated and easily 
identifiable. The MSTO is not an ideal distance indicator 
itself, since it exhibits considerable variation in absolute 
magnitude with differences in the age and metallicity of 
the stellar population, but it does provide a reasonable 
estimate of the stream's distance such that we can use it 
as an additional constraint on the position of the RC. 

Our distance measurements are thus made with a two- 
step process, whereby we first fit the bright edge of the 
MSTO, which unambiguously detects the Sgr stream but 
does not yield a precise distance, then use that measure- 
ment to set the range of apparent magnitudes which we 
use for detecting the red clump. The red clump fit then 
yields a precise distance. 

Towards this goal, we first present an overview of the 
PS1 survey in Section[2j and present the PS1 view of the 
southern Sgr stream, along with distance measurements, 
in Section [3l We discuss the results of those fits and 
their relation to other works in Section 21 and conclude 
in Section [5j 

2. OBSERVATIONS AND DATA PROCESSING 

The observations were conducted by the Panoramic 
Survey Telescope and Rapid Response System 1 (Pan- 
STARRS1, hereafter referred to as PS1). PS1 is a 1.8m 
telescope on Haleakela, Maui, which has been conduct- 
ing a multi-faceted survey since May 2010. The survey 
is designed both to detect transient, variable and moving 
objects, such as supernovae, Kuiper belt objects, aster- 
oids, and stellar transits, and also to provide data to a 
number of static-sky projects, such as large scale struc- 
ture, galaxy properties, and Milky Way structure. The 
survey covers the entire sky north of declination —30° 
(the so-called 3ir survey), with 12 specific fields targeted 
for more frequent and deeper observations (the medium- 
deep survey). 

The pattern of observations on the sky is organized 
by a fixed tessellation pattern which sets the boresight 
pointing locations. The timing of visits to each posi- 
tion is set by the need to obtain an adequate baseline 
in time for detection of transients and moving objects. 
The resulting set of observations for each position is de- 
signed to comprise of four observations of ^40 seconds 
in each of the five PS1 filters, per observing season. The 
images are processed by the Pan-ST ARRSl Image Pro- 
cessing Pipeline (IPP, [Ma gnierl [20061 ) . which handles all 
image processing steps including production of an object 
catalog for each image. Because of the computationally 
intensive nature of the processing on even single images, 
the IPP does not immediately produce stacked images. 
This is a necessary trade-off to enable the IPP to process 
images rapidly after the observations were taken, which 
is critical for the surveys of moving objects which require 
rapid follow-up observations. 

As a result of this processing strategy, the currently 
available 37r-survey catalogs contain only the results of 
photometry performed on single images. These catalogs 
are then cross-matched and merged to form an all-sky 



catalog. This does not increase the photometric depth 
in the same way that stacking images would, so it is not 
possible to detect objects fainter than the limiting mag- 
nitude of the deepest image. In order to reject spurious 
detections from instrumental artifacts we require that 
objects be detected multiple times, which makes our de- 
tection probability dependent on the number of visits to 
a position on the sky. 

The depth of the resulting catalog is therefore depen- 
dent on the number and depth of individual exposures in 
each band in a particular region of the sky. The number 
of visits to each position on the sky varies due to weather 
and telescope downtime affecting the scheduling of vis- 
its. There are also variations in observing conditions be- 
tween images, since the survey is designed to take data 
in marginal photometric conditions, which can affect the 
photometric completeness at faint magnitudes. To assess 
this variation, we have cross-matched the PS1 catalog, 
with the SDSS "stripe 82" coadd catalog ([Annis et al.l 
120111 ). The coadded SDSS catalog is better than 90% 
complete for stars down to 23rd magnitude in g, r, and 
i bands, so we can assume that the any non-detections 
of SDSS stars in the PS1 catalog is the result of PS1 
incompleteness. From this cross-matching we can estab- 
lish the 50% completeness magnitude over the range of 
observing conditions seen by the fields overlapping stripe 
82. In g-band the 50% completeness ranges from approx- 
imately gpi = 21.4 to 22.0, in r-band from rpi = 21.2 
to 21.8, and in i-band from i P1 = 21.0 to 21.8. While 
these depths are fainter the range of magnitudes we use 
for in this work, because the completeness is not a step 
function as a function of magnitude there will be some 
variation in completeness even in the brighter data we 
use for the MSTO fitting. From this cross-matching we 
can see that the completeness at ip\ = 21.0 varies by 
approximately 20% across the range of observing condi- 
tions in stripe 82. This variation is reduced to 5% or less 
at magnitudes of ipi = 19.0 or brighter. Much of this 
variation is likely to be related to our ability to separate 
stars from galaxies, which is progressively degraded un- 
der worse observing conditions. The photometric uncer- 
tainty on detected objects is less variable and frequently 
very small, on average ranging from 0.04 to 0.08 magni- 
tudes near the MSTO (at approximately 21.5 mag in g 
and r filters), and less than 0.03 mag at RC magnitudes 
(in r and i filters). 

One of the major challenges with any large survey 
is photometric calibration. The initial calibration of 
the IPP catalog is performed by referencing stars in 
each image to stars fr om the 2 Micron All Sky Survey 
(|Skrutskie et al.|[2006D . which provides a consistent way 
of calibrating single images across the entire sky. How- 
ever, now that the survey has obtained a sufficiently large 
set of observations with overlapping images, it is possi- 
ble to self-consistently calibrate the survey using stars 
observed in multiple overlapping observations. Com- 
monly called "iibercalib ration" , this calibration str ategy 
was adopted by SDSS (jPadmanabhan et alJ l2008f) an d 
has been implemented for PS1 by iSchlaflv et al.l (|2012|) . 
Comparisons between SDSS and the ubercalibrated PS1 
data show that the resulting measurem ents agree to 
< 10 mmag in the gpi, rpi, and ipifilters (|Schlaflv et al.1 
2012). This calibration is applied to the data using the 
Large Survey Database software (LSD, Juric, in prep.), 
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Figure 1. Map of color-selected Sgr MSTO stars, shown in galactic coordinates and centered o n the south galactic ca p. The right image 
shows the position of the bright arm (solid blue) and the faint arm (solid red) as measured bv lKoposov et all <2012f ). while the dashed 
versions show these lines extrapolated beyond the SDSS coverage. The position of the dwarf is indicated by the star on the right. The 
dashed rectangle indicates the region shown in Figure [2] 



which provides a fast and scalable interface to the large 
volume of data required. Analysis of the calibrated data 
was also performed using LSD. All of the observations 
were extinction corrected and de-redd ened using the red- 
dening maps of iSchlegel et al.l (|1998f) . We make use of 
all PS1 37T data obtained through January 18, 2012, at 
which time the survey has imaged the entire 3tt target 
area, though not always at the target photometric depth. 

3. THE PSl VIEW OF THE SGR STREAM 

Figure [1] shows a map of stars in the southern Galactic 
hemisphere with colors similar to that of the Sgr stream 
MSTO (0.0 < (gpi - rpi) < 0.4 and 20.25 < g Ph0 < 
22.10). The Galactic center is on the right edge, and 
the anticenter is on the left. The map is most sensitive 
to MSTO stars at heliocentric distances between 25 and 
45 kpc (assuming an old, metal poor population). The 
position of the bright arm of the stream is illustrated by 
the blue line on the right panel and the posit ion of the 
faint arm (as measured bv lKoposov et ai1l2012l ) is shown 
by the red line. The stream is clearly visible at b < 
—60°, nearly passing through the south Galactic pole. 
Towards the Galactic center, the stream skirts the edge of 
the PSl coverage and becomes lost in the higher density 
environment of the Galactic center. In the anticenter 
direction, a combination of varying observation depth 
and the increasing heliocentric distance of the stream 
makes it difficult to discern the stream at b > —45°. 

The Sgr stream in the s outhern hemisphere has also 
been seen in SDSS data (IKoposov et all 12011) over a 
somewhat smaller area. The SDSS and PSl coverages 
are shown in Figure O IKoposov et al.l (|2012h have also 
shown that the Sgr stream in the southern Galactic hemi- 
sphere is bifurcated in much the same way as the north- 
ern part of the stream. They show that the density of 
MSTO stars is asymmetric, with MSTO stars extending 
further perpendicular to the stream on the northern side 



(top of Figure[2]) than the southern side. This is apparent 
in both the SDSS maps and the PSl maps. 

3.1. Distance Measurements 

We divide the length of the stream visible in the MSTO 
map into ni ne regions space d equa lly along the Sgr plane 
defined by iMajewski et al.l (|2003l ) . Though this plane 
was designed to fit the position of the stream on the 
sky, the peak of the MSTO star density we observe lies 
slightly to the south of this plane. In order to improve 
the signal to noise of our measurements, for each of the 
nine regions we construct a histogram of MSTO star den- 
sity as a function of distance off of the Sgr plane, then 
center our target on-source regions on the peak of these 
histograms. In general the region centers are approxi- 
mately one degree south of the Sgr plane, but each region 
is fit individually so there is some variation. Not all of 
this variation is necessarily physical though, since we are 
also affected by small gaps in the PSl sky coverage. As 
a result, the positions of these regions may only approxi- 
mate the true position of the Sgr stream on the sky. The 
on-source regions cover 6° along the stream and are 5° 
wide perpendicular to the stream. The field centers are 
listed in Table [1] and the on-source regions are shown in 
the top panel of Figure [2] in blue. 

In addition to placing on-source regions near the Sgr 
plane, we also placed regions 17° away from the plane 
to serve as background regions for the MSTO fits. Ide- 
ally these background regions would be located at similar 
galactic latitude to the target regions, but at high galac- 
tic latitude the PSl sky coverage does not include suffi- 
cient area that is not contaminated with the Sgr stream 
itself for this to be possible. Instead we adopted the 
strategy of background regions parallel to the stream, 
which for regions near the disk tends to approximate a 
constant galactic latitude selection process due to the po- 
sition of the stream on the sky. At high galactic latitude 
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Figure 2. Zoomed-in view of Figurc[T] focusing on the Sgr stream. 
The top panel shows the data in PS1, while the bottom shows the 
same region in the SDSS. The top panel also shows the regions 
targeting the bright arm of the stream in blue, the faint arm regions 
in green, and the background regions in red. The coordinate grid 
labels show Galactic longitude, and the concentric rings are each 
separated by 15° Galactic latitude. 

this approximation begins to break down, but the back- 
ground contamination is much less significant since those 
fields are far away from the disk. We also placed regions 
that targeted the faint arm of the stream approximately 
10° away from the plane, where the MSTO distribution 
shows a second peak or at least show signs of being spa- 
tially extended to the north. The resulting regions are 
shown in the top panel of Figure [2] with the faint arm 
regions in green and the background regions in red. 

A background-subtracted color-magnitude diagram 
(CMD) for a sample region is shown in Figure [3] The 
background region was scaled to match the area of the 
target region. An isochrone of an 11 Gyr old po pu- 
lation with [Fc/H]= -1.2 from iMarigo et all (|2008l ) is 
shown for ref erence (computed with filter curves for the 
PS1 system; iTonrv et ail 1201 2f > . This metallicity and 
age match measurements of the stream in the north, 
thou gh the stream ex hibits a broad range of metallici- 
ties ()Chou et al.ll2007l) . The Hess diagrams for the target 
fields were used to determine the color cuts that would 
best select the RC (0.17 < (r P1 - i P1 ) < 0.23) and the 
MSTO (0.0 < (rpi - ipi) < 0.10). The RC color selec- 
tion is shown by tall white rectangle. 

These color cuts were then used to construct his- 
tograms of the apparent magnitudes of stars with the 
selected colors in each region. These regions are shown 
in Figure [4] for the bright arm regions and Figure [5] 
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Figure 3. CMDs for region E, showing the original on-source 
CMD (loft), the background-subtracted CMD (center), and the 
background-subtracted CMD (right) with the RC color-cuts and 
an isochro ne for an 1 1 Gyr old population with [Fc/H]= —1.2 over- 
laid ffrom lMarigo et al.ll200§h . The horizontal arrow indicates the 
apparent magnitude of the red clump. 

for the faint arm. For presentation purposes the fig- 
ur es show "running his togra ms", similar to t hose used 
in iCorrenti et all (j2010f ) and iBellazzini et all ()2005h . In 
these running histograms, each point is spaced by 0.02 
magnitudes (the step size), but the value of each point is 
the number of stars within a 0.2 magnitude bin centered 
on the point (the bin size). The procedure is essentially 
a boxcar filter. The resulting step values are not statis- 
tically independent of each other, but the running his- 
togram has the advantage that the visual appearance of 
the data is not affected by the particular bin positions 
in the histogram. That is, if the RC happens to fall on 
the edge between two bins of a conventional histogram, 
its signal appears more spread out than if it fell in the 
center of a single bin. The bin width is also chosen so 
that it roughly corresponds to the size of the RC, thus 
maximizing the signal to noise. The running histograms 
are primarily a presentation tool; all of the statistical fits 
which made use of histograms used "normal" , fixed-bin 
histograms. 

3.1.1. MSTO Fit 

After constructing the histograms of the MSTO and 
RC stars, we fit the apparent magnitude of the MSTO 
using a maximum likelihood method, modeling the edge 
of the MSTO as a step function superimposed on smooth 
background, as measured from the off-stream field and fit 
by a third-order polynomial. The step function was then 
convolved with a Gaussian with a width of a = 0.10 
magnitudes to account for both photometric scatter and 
intrinsic scatter in the stellar luminosities. The result of 
these fits can be seen in the insets of Figure HI In general 
these fits are very satisfactory, but it is clear that the 
edge of the MSTO is not a very sharp transition; the 
transition can span as much as half of a magnitude. 

To use the MSTO fit as a prior on the RC fit, we 
must determine the difference in apparent magnitude be- 
tween the two features. This calibration is inexact due 
to the various complicating factors associated with the 
MSTO but still yields useful results, since for this ap- 
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Figure 4. Running histograms for the regions on the bright arm of the stream. The blue line shows the running histogram, while the 
black lines show the best fitting model (or multiple models in the case of multiple candidate peaks.) The horizontal dashed line indicates 
the level at which an overdensity would have a 3<r significance over the computed background level. The insets shows the MSTO fits used 
as a prior on the RC fits, with the on-stream data shown in blue, an off-stream background region shown in red, and the best fit model in 
black. The arrow shows the position of the red clump suggested by the MSTO fit, and the range of RC values allowed by the MSTO prior 
is shown by the vertical dotted lines. In cases where there was no peak with significance greater than 2-cr, we do not plot a fit and instead 
denote the background level with a solid horizontal line. 

A necessary step in fitting the RC is the determination 
of the smooth background level, which by number dom- 
inates over the RC. As shown in the CMD (Figure [3]), 
the RC color cut we used is just slightly redder than 
the region of the CMD populated by nearby halo stars, 
which significantly outnumber the RC stars we are at- 
tempting to detect. Though the peak density of these 
stars occurs slightly bluer in color than the RC, both 
intrinsic scatter in color and photometric error causes 
some of these stars to fall into our color cuts. The ap- 
parent magnitude distribution of stars in this region of 
the CMD is approximately constant and not strongly af- 
fected by observing conditions. As a result, we model 
the RC background as a third order polynomial in color 
(between (rpi — ipi)o = 0.15 and 0.35) and a constant 
in apparent magnitude. The polynomial fitting excludes 
the RC color window to prevent biasing the background 
estimation. 

This approach to modeling the background has a num- 
ber of advantages. The primary reason for adopting this 
method is that the amount of background contamination 
to the RC color cut is highly sensitive to photometric 
errors in the observations. Greater observational uncer- 
tainties causes the large number of MSTO-color stars to 



plication we are only using it for different parts of the 
same stream. Since the composition of stellar popula- 
tions in the stream is unlikely to vary greatly across the 
relatively modest length of stream available to us, the 
RC-MSTO offset should also only exhibit small varia- 
tions. We com puted this off s et by creating synthetic 
C MDs from the | Marigo et all (|2008l ) isochrones (using 
a iChabrieil (|200ll ) initial mass function), then apply- 
ing the same color cuts and using the same maximum 
likelihood fits for the MSTO and RC position as were 
used on the PS1 data. Using an isochrone for an 11 
Gyr old, [Fe/H] = -1.2 population the RC-MSTO offset 
was computed to be 2.50 magnitudes in i-band. Varying 
the metallicity of the isochrone from [Fe/H] = —0.5 to 
[Fe/H]= —1.3 (roughly the range of stream m etallicitics 
seen in the northern arm bv lChou etafl (|2007l )) and also 
varying the age of isochrone from 6 Gyr to 11 Gyr, pro- 
duced RC-MSTO offsets ranging from 2.71 through 2.40. 
The prior we adopt on the RC fit is therefore broader 
than this range of variation, but centered on an offset of 
2.50 magnitudes. 

3.1.2. RC Fit 
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Figure 5. RC fits for the regions placed on the faint arm of the stream. The lines are the same as in Figure [4] However, because we set 
the window for the bifurcation fits with the measured RC position in the on-stream regions, the allowed RC range (vertical dotted lines) 
may not align exactly with the center position suggested by the MSTO fit (vertical arrow). 



occupy a wider range of colors, causing more to fall into 
the RC color cut and hence a higher background. This 
is particularly true for PS1, which has data spanning a 
range of observing conditions. This makes it very dif- 
ficult to use off-stream regions as representative of the 
background level, since it is practically impossible to find 
suitable background regions that reproduce the observ- 
ing conditions of an on-stream region. For this reason we 
also do not attempt to fit a precise function to the back- 
ground as a function of apparent magnitude; though in 
several of the regions the background deviates from flat- 
ness, it generally does not do so in a way that could be 
modeled by comparison to an off-source region. Addi- 
tionally, since our RC fits are constrained by the MSTO 
prior to a small range in apparent magnitude, it only 
matters that our background estimation matches that 
small part of the histogram well. The assumption of a 
constant background level is quite suitable under these 
conditions. 

After we have obtained the prior on the RC fit from 
the MSTO, we fit to the data the sum of a Gaussian 
function with a fixed width of 0.1 magnitudes and the 
constant background level. The MSTO fit and the RC- 
MSTO calibration were used to set a flat prior within 0.6 
mag of the center. The prior is broad enough to account 
for any variation in the age and mctallicity of the stream, 
along with the inherent uncertainty and variability in the 



MSTO measurement as a result of observational factors. 
The scale factor on the Gaussian is then marginalized 
out, and the best fitting magnitude and uncertainty is 
found by taking moments of the likelihood. In general 
this produces good fits, with the statistical uncertainty 
in the measurement significantly smaller than the uncer- 
tainty of the RC calibration. 

In a small number of fields, multiple distinct peaks ap- 
pear in the histograms with significance greater than 3c 
In many of these cases one of the peaks is more likely 
given the position of the stream in adjacent fields, but 
for completeness we report both distances. In some other 
fields, the RC appears to span a larger range of abso- 
lute magnitudes than expected, or exhibits some asym- 
metry that may indicate some feature is present other 
than the "pure" Sgr RC at a single distance. In these 
cases we do not have enough information to report mul- 
tiple distances, since it would be an ambiguous exercise 
to decompose the noisy histogram into several RC com- 
ponents. Instead we report a single distance but with 
an uncertainty that reflects the entire range of magni- 
tudes where the value of the running histogram exceeds 
the background by 3c, as estimated from pure Poisson 
noise (corresponding to the horizontal dashed lines in 
Figures @] and [3]) . This effectively determines the range 
of magnitudes at which, if our assumption of a constant 
background holds, a statistically significant overdensity 
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exists, even if we are unable to conclusively identify its 
origin. 

An additional source of uncertainty comes from our 
assumption of the faint arm's position on the sky. In 
selecting our regions for the distance measurements we 
have assumed that the faint arm parallels the bright arm 
over the region seen by PSl, but our MSTO maps are 
unable to localize the faint arm pr ecisely enough to verif y 
that this is the case. As shown bv lKoposov et al.l (|2012fl . 
the streams may be converging towards the Sgr dwarf 
at a rate of 0.05 degrees per degree along the stream. 
Our target regions are designed to be sufficiently wide 
perpendicular to the stream to encompass such gradients, 
but we have also tested our faint arm fits with the target 
regions placed along paths converging or diverging from 
the bright arm, with angles ranging from 0.10 deg/dcg 
to -0.10 deg/deg. The results still generally agree well 
with the results from the parallel-stream assumption, but 
there is some variation in distances of approximately the 
same order or less than the uncertainties on individual 
distance measurements. As a result we have incorporated 
these additional uncertainties in our reported errors. 

In order to use the RC as a distance indicator it is 
necessary to have a calibration of the feature's absolute 
magnitude. This is complicated in the case of the RC 
by the fact that its absolute magnitude is weakly de- 
pendent on the age and metallicity of the ste llar pop- 
ulation (JsHdeTetj^ [ljjH [jimene^iFanillll). These 
variations can be minimized by using a red photometric 
passband (jPaczvnski fc Staneklll998l f and by calibrating 
the RC magnitude with stars of the same stellar popu- 
lation at a known dista nce. This "differentia l" measure- 
ment was employed by iCorrenti et al.l (|2010l ) and is the 
method we use here. Since the Sgr dwarf lies within the 
PSl survey area we are able to directly use the survey's 
observations of the main body of the dwarf to calibrate 
our RC distance measurements along the stream. This 
avoids the need for any conversion between photomet- 
ric systems. The histogram of RC color-selected stars 
in the main body of the dwarf is shown in Figure |6l in 
which the RC is immediately obvious. We fit the Sgr 
dwarf RC using the same methods as for the stream, but 
with a polynomial background measured from a neigh- 
boring off-source region instead of the constant used for 
the stream. The resulting fit measures the apparent mag- 
nitude of the RC to be i P1 = 17.27 ± 0. 05. We adopt 
the Sg r dwarf distance measurement by I Monaco et alJ 
(|2004f ). which used the tip of the red giant branch to mea- 
sure a distance of 26.30 ± 1.8 kpc, o r a distance modulus 
of (m -M) a = 17.10±0.15 (but see lKunder fc Chaboverl 
(|2009f ) for a recent alternative distance). All of our dis- 
tance measurements are based on this number, and so 
all of our reported distances will scale directly with any 
revised value of the distance to the dwarf. 

4. RESULTS & DISCUSSION 

Our measured distances to the stream are shown in 
Figure [7j with the bright arm shown in red points and 
the faint arm in blue. While the faint arm is just on 
the edge of detectability in several individual fields, the 
consistent behavior of the faint arm across the length 
of the stream does suggest that these detections do 
form a coherent structure. Despite the significant un- 
certainties on each individual measurement of the dis- 
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Figure 6. Running histogram of RC color-selected stars in the 
main body of the Sgr dwarf. The solid line shows the histogram 
from the dwarf, while the dashed line shows the histogram for a 
background region. 

tance, the faint arm appears 3-5 kpc closer to the sun 
than the bright arm in several different fields. A dis- 
tance offset between the two arms of the stream is also 
seen in the northe r n Galactic hemisphere, as see n by 
Yannv et al.l (l2009l).lNiederste-Ostholt et al.l (|2010h . and 
Ruhland et alJ (|2011D . where the faint arm (also called 
the "B" stream) also lies slightly closer to the sun than 
the bright stream (the "A" stream). 

The distance to the main body of the trailing arm of 
the Sgr stream has also been measured with the red 
clump in SDSS data by IKoposov et alJ (|2012l) . Their 
distances are shown in Figure [7] as the solid line. Com- 
paring their results to ours for the bright arm of the 
stream, there is a clear difference of 7-10 kpc, which is 
consistent along the length of the stream. This appears 
to match our measured distances t o the faint arm, but we 
believe that this is a coincidence. IKoposov et alJ ([2012D 
show the density of RC stars along the stream in SDSS 
as a function of apparent magnitude (their Figure 5), 
and examination of this reveals that their detections of 
the RC lie at nearly the exact same apparent magni- 
tudes as ours. However, IKoposov et al.l (|2012l ) adopt 
an absolute ma g nitud e of the RC as Mi = 0.6, citing 
iBellazzini et alJ (|2006| ). while we use M i:PS = 0.17 (the 
SDSS and PSl i-band filters are nearly identical). This 
discrepancy is significantly larger than the uncertainty on 
the distance to the Sgr dwarf, which is 0.15 magnitudes. 
In Figur e [71 the dashed l i ne sho ws the distance measure- 
ments of IKoposov et all (|2012f) shifted by the difference 
between these calibrations, and the resulting distances 
agree well with our data. It thus appears that the dis- 
crepancy between the two results is caused by differences 
in the two RC absolute magnit ude calib rations. 

Tho ugh lKoposov et al.l ([2012D refers to lBellazzini et ail 
(|200l as the source of their i-band RC calibration, it is 
unclear how this i-band measurement was obtained, since 
IBellazzini et~aTl (|2006| ) uses observations performed en- 
tirely in the B and V bands and include no mention of 
the SDSS-i filt e r. W ith no further details provided by 
IKoposov et alJ (|2012D we cannot diagnose the process 
that lead to this M* = 0.6 calibration. We are confi- 
dent that our own calibration is robust: we detect the 
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RC in the Sgr dwarf itself, using identical filters and ob- 
served as part of the same survey as our detections of the 
stream, thus minimizing the systematic uncertainty as- 
sociated with color transformations and comparisons be- 
tween data of different origins. One other possible source 
of variation between calibrations is the assumed extinc- 
tion values to the Sgr dwarf, since it lies close to the 
Galactic plane, but it is unlikely that extinction could 
account for the magni tude of the calibration difference. 
iBellazzini et al.1 (|2006l ) computes an average reddening 
of E(B-V) =0.1 16 , while we use the value from the maps 
of lSchlegel et all (|1998f ) which report E(B-V)=0.15. In 
terms of extinction this amounts to a difference in Ai of 
0.07 magnitudes, much less than the difference in the two 
calibrations. 

In Figure El the measured distances are compared to 
the numerical mo d el of t he Sgr tidal streams produced by 
lLaw fc Majewskil (|2010t l. The fi gure shows a projection 
of the stream onto the Galactic X-Z plane, which approx- 
imates the plane of the Sgr stream, and the Y-Z plane. 
The agreement between our distances and the position of 
the trailing (southern) arm of the Sgr debris is excellent. 
This is in light of the fact that the simulation was not 
selected to match the distance to the stream, but only 
the measured radial velocities to the leading and trailing 
streams. In the right panel of Figure [51 the Sgr dwarf 
appears slightly offset from the plane of the bright arm 
of the stream, but this is the result of a projection effect. 
Because the Sgr stream does not lie exactly in the X-Z 
plane, the projection of the stream onto the Y-Z plane 
will exhibit some variation in the Y direction. Because 
the dwarf lies on the opposite side of the Galactic center 
from the majority of the stream points, they will appear 
to be offset in Y when projected in this manner. 

Figure|8jal so shows the pos i tion o f the Cetus stream, as 
measured bv iNewberg et al.l (120091 ) . Though the stream 
crosses the position of the Sgr stream on the sky, and 
in the X-Z plane of Figure [5J (left panel) appears to be 
coincident with detections of the faint arm, the X-Z pro- 
jection of Figure [8J clearly shows that the Cetus stream 
passes behind the Sgr stream and not through it. It is 
therefore unlikely that our detections are confused with 
the Cetus structure. 

Though there is currently no model for the origin of 
the northern and southern bifurcations that matches all 
of the available data, it is still instructive to examine 
how our results compare to the existing models in order 
to develop a sense of how a plausible bifurcation may 
behave in simulations. One hypothesis is that the Sgr 
dwarf was originally a disc galaxy, and that the material 
in the faint arm was stripped during a different peri- 
centr ic passage than the bright arm (jPenarrubia et ahl 
120101 ). This theory is now thought to be untenable, since 
the S gr dwarf is observed to have very little residual rota- 
tion (IPeharrubia et al.ll2011l ). but the simulations used in 
iPeharrubia et al.l ( 20101) are still useful for assessing the 
difficulty of reproducing the southern bifurcation simul- 
taneously with the north. In the scenario that best repro- 
duced the northern bifurcation, where the two streams 
are nearly coincident, the material from the two streams 
in the southern hemisphere is widely divergent. Material 
stripped in the same pericentric passage as the north- 
ern faint arm rapidly increases in helicentric distance to- 
wards the Galactic anticenter (from roughly 25kpc to 
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Figure 7. Heliocentric distances to the bright arm (red, filled cir- 
cles) and faint arm re gions (blue, open circle s), plotted in the Sgr 
coordinate system of[Maicwski et al. (2003). The smaller points 
indicate detections with significances between 2 and 3cr. The 
solid line indicates t he distances to the Sgr stream reported by 
IKoposov et a l. (2012), w hile t he dashed line shows the measure- 
ments of Koposov et al. (2012) but adjusted to our calibration of 
the RC absolute magnitude (Mi = 0.17). The Sgr dwarf is located 
at As 9 r,o = 0. 

80kpc) , while material from the bright arm stays at rela- 
tively constant heliocentric distance (from approximately 
20kpc to lOkpc). These simulations were not designed to 
reproduce the stream's behavior in the southern hemi- 
sphere and cannot be faulted for not doing so, but the 
discrepancy highlights the difficulty inherent in the prob- 
lem. 

5. CONCLUSIONS 

In this work we have used data from Pan-STARRSl to 
show the spatial extent of the Sgr stream over 60° of its 
orbit in the southern Galactic hemisphere. The position 
we observe of the stream on the sky matches well wit h ob- 
servations from o ther surveys (|Koposov et al.ll2012T ) and 
with simulations (|Law &: Majewski 2 0101 ). We have also 
shown that the stream in the southern hemisphere ex- 
hibits a bifurcation similar to t hat seen in the northe rn 
hemisphere, as was reported bv IKoposov et all (|2012| ). 

Using a combination of the MSTO and the RC as dis- 
tance indicators, we have measured the distance to both 
arms of the southern stream. Our resu lts for the stream 
again agree well with the simulations of lLaw fc Maiewskl 
(I2010D. but disagree with the distances measured by 
IKoposov et al.l (120121 ) . We believe that this disagreement 
is the result of differences in calibrations of the RC ab- 
solute magnitude, and not differences between the SDSS 
and PS1 surveys, but it is not possible for us to further 
deduce the exact cause of the disagreement. Our cali- 
bration strategy, based on direct observations of the Sgr 
dwarf within the PS1 survey, gives us confidence that we 
have one of the most direct calibrations possible, free of 
any photometric transforms between data sources. Our 
distance measurements consistently place faint arm of 
the stream slightly closer in heliocentric distance than 
the bright arm. This is similar to the behavior seen in 
the northern hemisphere, again suggesting that the bi- 
furcation is the result of some intrinsic behavior in the 
accreting Sgr system and not a coincidental overlap of 
multiple wraps. 

Unfortunately there is no model for the bifurca- 
tion that can adequately explain the observed behav- 
ior. Models requiring multiple wraps of the stream 




Figure 8. Detected positions of the Sgr stream projected onto the Galactic X-Z plane (left), which approximates the orbital plane of the 
stream, and in the Y-Z plane (right), which is perpendicular to the orbit. The position of the Galactic disk is shown by the horizontal 
line, and the sun is marked with an X. In the south, red circles and blue triangles signs indicate the position of statistically significant 
detections of the RC on the bright arm and the faint arm fields, respectively. The s maller points indicate d e tectio ns with significances 
between 2 and 3cr. In the northern hemisphere, SDSS detections of the stream from Nicdcrstc-Ostholt ct al. (2010) arc shown with red 
circles for the bright arm, and blue triangles again for the faint arm. The position of the Cetus stream UNewberg et al. 2009) is also shown 
by the green squares. The black dots are data from the simulations of Law & Majcwski (2010). 



Field 

Field center (RA.Dec) 
Field center (A s3riQ ,B sgr ^) 
Distance (kpc) 


A 

(50,8) 
(126,-1) 
35.6 


B 

(43,5) 
(118,-1) 
34.2,27.8 


C 
(37,1) 
(110,-1) 
33.1,29.5* 


D 

(29,-2) 
(102,0) 
29.5 


E 

(23,-7) 
(94,-0) 
27.4 


F 

(16,-11) 
(87,-1) 
26.7 


G 
(9,-14) 
(79,-1) 
26.7* 


H 

(3,-18) 
(72,-1) 


I 

(356,-20) 
(65,-0) 


Bifurcation center (RA,Dcc) 
Bifurcation center (A sgr: Q,B Bgry Q) 
Bifurcation Distance (kpc) 


(47,18) 
(128,9) 
42.5*,27.8*,32.7* 


(40,14) 
(120,9) 
28.2,23.0 


(32,11) 
(112,10) 
26.3,31.8* 


(25,7) 
(103,10) 
23.3* 


(18,3) 
(95,10) 
26.3* 


(11,-1) 
(87,10) 


(4,-5) 
(79,10) 
24.1* 


(358,-8) 
(72,10) 
33.6* 


(351,-11) 
(64,10) 



Table 1 

Positions of the regions along the Sgr stream used to detect the RC, along with the corresponding distance measurements. In cases where 
multiple overdensities could correspond to the red clump, the distances from each candidate position are listed. Detections with 

signifiancc between 2 and 3 — cr arc marked with an asterisk. 



(jFellhauer et al.l 120061 ) or internal rotation of the Sgr 
dwarf appeared to be untenable even prior to the discov- 
ery of the bifurcation in the south (jYanny et al.1 120091 : 
iNiederste-Qstholt et al.ll2010 ). It h as been speculatively 
suggested bv lKoposov et a!T ( 2012D that what we call the 
Sgr stream could be the result of the accretion of a pair of 
dwarf galaxies simultaneously. Our observations should 
provide a quantitative basis against which simulations 
of a double accretion scenario can be tested. It is re- 
markable that one of the most prominent and most well- 
studied tidal streams has proven to be the most elusive 
streams to explain. While our observations can at the 
moment only deepen the mystery, we are hopeful that 
the added information on the stream's behavior will be 
used to properly explain the bifurcation in future efforts. 
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